Full quantum reconstruction of vortex states 
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We propose a complete tomographic reconstruction of vortex states carrying orbital angular momentum. The 
scheme determines the angular probability distribution of the state at different times under free evolution. To 
represent the quantum state we introduce a bona fide Wigner function defined on the discrete cylinder, which is 
the natural phase space for the pair angle-angular momentum. The feasibility of the proposal is addressed. 
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Phase singularities of wave fields were brought to atten- 
tion in a number of seminal papers HI 0. Ht] . They manifest 
in electron wave packets 10], quantum Hall fluids J5h, super- 
media fl, ferromagnets U, Bose-Einstein condensates JH, 
acoustical waves Jgfl, and light fields Jioll . to cite some rele- 
vant examples. Special attention has been paid to the particu- 
lar case of vortices, which correspond to helicoidal waves ex- 
hibiting a pure screw phase dislocation along the propagation 
axis; i.e., an azimuthal phase dependence e ^. This means 
that I plays the role of a topological charge: the phase changes 
its value in £ cycles of 2ir in any closed circuit about this axis, 
while the amplitude is zero there. 

One of the most interesting properties of vortices is that 
they carry orbital angular momentum (OAM): indeed, the in- 
teger £ can be seen as the eigenvalue of the OAM operator and 
its sign defines the helicity or direction of rotation. In fact, 
OAM can be easily transferred, as demonstrated in a number 
of recent challenging experiments with optically trapped mi- 
croparticles 111 111 . 

At the microscopic level, quantized vortices are now rou- 
tinely produced by, e.g., mechanical stirring of ultracold 
atomic clouds. Remarkable images of large lattices containing 
hundreds of vortices in an Abrikosov-type triangular configu- 
ration have been obtained fl2| . However, for a complete infor- 
mation about the corresponding state, one needs to go beyond 
this mere photographic picture and perform a full tomography. 
Efficient methods of state reconstruction are of the greatest 
relevance for quantum optics. Since the first theoretical pro- 
posals, this discipline has witnessed significant growth 11311 
and laboratory demonstrations of state tomography are numer- 
ous and span a broad range of physical systems. The essence 
of vortices is their helicoidal structure, while other aspects, 
such as energy distribution, may be often ignored for many 
purposes: in this sense, a fully quantum reconstruction of the 
vortex content is pending. 

Any reliable quantum tomographical scheme requires three 
key ingredients 11411 : the availability of a tomographically 
complete measurement, a suitable representation of the quan- 
tum state, and a robust algorithm to invert the experimental 
data. When these conditions are not met, reconstruction be- 
comes difficult, if not impossible: this has been the case so 
far for vortex states. The construction of a proper Wigner 



function for them (or any other quasiprobability distribution 
in phase space) is still an open question. Although some in- 
teresting attempts have been published Il5l llfjll . they seem 
of difficult application for the problem at hand. The twofold 
goal of this paper is precisely to fill this long overdue gap. We 
will provide a simple Wigner function with a clear geomet- 
rical meaning and also a universal tomographic reconstruc- 
tion scheme. One could think in using the standard Wigner 
function for the transverse coordinate-momentum variables 
instead. However, this is not the right answer in quantum 
mechanics: such a representation conveys redundant infor- 
mation, it is impossible to plot, and it hides all the relevant 
angular information. 

We consider rotations by an angle <fi generated by the OAM 
operator along the z axis, which for simplicity we shall de- 
note henceforth as L. We do not want to enter in a long and 
sterile discussion about the possible existence of an angle op- 
erator. For our purposes here the simplest choice is to use the 
complex exponential of the angle E = e -4 *, which satisfies 
the commutation relation [E, L] = E. The action of E on 
the OAM eigenstates is E\£) = \£ — 1), and it possesses then 
a simple implementation by means of a phase mask remov- 
ing a unit charge lfl7ll . Since the integer I runs from — oo to 
+oo, E is a unitary operator whose normalized eigenvectors 
|0) describe states with well-defined angular position. In the 
representation generated by them, L acts as — id^ (in units 
h = 1). Note in passing that one could intuitively expect a 
Fourier-like relationship between angle and OAM, which can 
be expressed in this context as e~ 1 ^ L \4>) = \<f> — <f>'). 

For the standard harmonic oscillator, the phase space is the 
plane K 2 . In the case of angle and OAM the phase space is the 
discrete cylinder S 1 x Z (where S 1 is the unit circle). While 
quasiprobability distributions on the sphere are commonplace 
in quantum optics, to the best of our knowledge their counter- 
parts on the cylinder have never been used to describe angular 
variables. Given the key role played by the displacement op- 
erators in defining the Wigner function in the plane, we intro- 
duce a unitary displacement operator L>(£, (f>) on the discrete 
cylinder as 



(1) 
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where a(£,4>) is an undefined phase factor. Apart from 2ir- 
periodicity in <fi, the requirement of unitarity imposes the con- 
dition a(£, 4>) + a(— £,—(/)) = —£<f>. The displacement op- 
erators form a non-Hermitian orthogonal basis on the Hilbert 
space, in the sense that 



Tr[D(£, <j>) D\e, </>')] = 2ttSw S 2v {<t> - (/>') 



(2) 



where 82^ represents the periodic delta function (or Dirac 
comb) of period 2n. 

Next, we introduce the Wigner kernel as a kind of double 
Fourier transform of the displacement operators 



w(£,<t>) 
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d(j)' exp[-i(e'(j) - £<j>')]D(£', <j>') 



(3) 

where the integral extends to the 2ir interval within which the 
angle is defined. The integral of w(£, (j>) over the whole phase 
space yields unity and we may thus regard w(£, (f>) as equiva- 
lent to the phase-point operators introduced by Wootters 1 1811 . 
In addition, one can check that 



w(£, <j>) = b(i, <j>) w(o, o) b\i, 4>) . 



(4) 



Since for the plane and the sphere, the Wigner kernel can be 
seen as the transform of the parity by the displacement op- 
erators (3, this seems to call for interpreting w(0, 0) as the 
parity over our phase space. 

We next define the Wigner function of a quantum state de- 
scribed by the density matrix g as 



W{£,(t>) = Ti[Qw{l,4>)\, 



(5) 



Using the previous results for w{l, (j>) one can show that 
W(£ 1 4>) fulfills all the properties required for a reasonable 
interpretation as a quasiprobability distribution. Indeed, due 
to the hermiticity of the Wigner kernel, W(£, cj>) is real. It 
also provides the proper marginal distributions and it is co- 
variant, which means that if the state g' is obtained from g 
by a displacement in phase space g' = D(£' ', </>') g D'(£', </>'), 
then the Wigner function follows along rigidly: W'(£, cf>) = 
W(£ — <j) — <j>')- All these properties are fulfilled indepen- 
dently of the choice of the phase a(£, <fi). 

Since the displacement operators constitute a basis, we can 
write the expansion 



(6) 



= E / dtj)Q{l,<j>)i){l,<t)) 



whereg(^,0) = Tr[£ 0)]/(2tt). In terms of </>), the 
Wigner function has the representation 

= ^ E / W> . (7) 



To work out explicit examples, one needs to fix once for all 
the function a(£, cf>). One natural option is to set a(£, </>) = 



and then the Wigner kernel © reduces to 

e- m '*\l + £')(£- 1'\ 
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x e- li - 2t ' +1 ^\£" + £' + !){£" -£'\ 
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which agrees with the kernel derived by Plebahski and 
coworkers |2~ofl in the context of deformation quantization on 
the cylinder. 

For an OAM eigenstate \£q), we obtain 



(9) 



which is a quite reasonable Wigner function: it is flat in (f> and 
the integral over the whole phase space equals unity, reflecting 
the normalization of \£o). 

For an angle eigenstate \4>o), we get 

wW*,0) = -^<M0-0o). do) 

Now, it is flat in the conjugate variable I, and thus, the inte- 
gral over the whole phase space is not finite, which is a con- 
sequence of the fact that the state |0o) is unnormalized. 

The coherent states |£o, 4>o) (parametrized by points on the 
cylinder) introduced in Ref. 112 ill (see also Ref. 12211 for a de- 
tailed discussion of the properties of these relevant states) sat- 
isfy 



(£\£ ,<f>o) = 
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where #3 denotes the third Jacobi theta function. One imme- 
diately finds that 

e a(0/2-0 o ) 



which gives 
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FIG. 1 : Plot of the Wigner function for a coherent state with £o = 
and <j)o = 0. The cylinder extends vertically from t = — 4 to I — +4. 
We show the two corresponding marginal distributions. 



In spite of the fact that these states may be regarded as a basic 
set to construct the Wigner kernel, their Wigner function has 
no simple closed form, but it splits into two different functions 
with periods tt and 2tt, respectively. 

In Fig. 1 we show the Wigner function for the coherent state 
\£o = 0, 4>o = 0) plotted on the discrete cylinder. We can see 
a pronounced peak at <p = for £ = and slightly smaller 
ones for I = ±1. A closer look at the picture reveals also a 
remarkable fact: for values close to <f> = ±ir and £ = ±1, the 
Wigner function takes negative values. Actually, a numeric 
analysis suggests the existence of negativities close to <j> — 
±7r for odd values of I. We also plot the marginals obtained 
from Eq. ( fl~3l by integrating over (j> or summing over t. 

As our last example, we look at the superposition states 
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which have been proposed for applications in quantum exper- 
iments 12311 ■ The analysis can be carried out for more general 
superpositions, but (fl3T l is enough to display the relevant fea- 
tures. The final result is 



1 

47T 



(Su +8e-i ) 



1 

'2tt 



cos(</>o-8c/>) S eo ■ (16) 



The Wigner function presents then three contributions: two 
flat slices coming from the states \£ ) and | — £ ) and their 
interference located at the origin (see Fig. 2). 

To complete our theory we also propose a reconstruction 
scheme for these sates. A reconstruction of g (or, equiva- 
lently, of its Wigner function) is tantamount to finding the 
coefficients g(£, cf>). To this end, we need a tomographical 
measurement that allows us to reconstruct g(£, 4>). For £ = 0, 



FIG. 2: Plot of the Wigner function (and the corresponding 
marginals) for a superposition state H51 with £o = 3 and (f>o — n. 
and its marginal distributions. The cylinder extends vertically from 
£ = -4 to £ = +4. 



the coefficients g(0, (f>) read as 



(17) 



where we have made use of the fact that the undetermined 
phase e~ ta< - '^ can be set to 1 for all values of 4>, for it corre- 
sponds to displacements along one of the coordinate axes and 
no additional phase should be acquired. 

For £ 7^ 0, we introduce the operator U t = exp(— itL 2 /2). 
If we recall that the generic Hamiltonian of a quantum ro- 
tator is H = L 2 /2I (I being the rotational moment of in- 
ertia), the operator U t is simply the free propagator of the 
system (in appropriate units). This quantum rotator Hamil- 
tonian describes a variety of situations, such as molecular ro- 
tations i24ll . single-photon OAM |25] or the azimuthal evolu- 
tion of optical beams 12611 . For all of them, our proposal for 
u}{4>, t) is precisely to measure the angular distribution after a 
free evolution t, that is, 



u{<P,t) = m)\gm)) = QUt 



(18) 



In other cases, the scheme also works appropriately provided 
one can experimentally implement the action of L 2 . Using the 
representations © for g and (Q~|) for D, g(£, tp) turns out to be 

g{£, </>) = -L e --(^) f dtf e- u * w^, <j>/i) . (19) 

27T J27T 

Plugging these coefficients into Eq. (0, we get the reconstruc- 
tion of the Wigner function as 



(20) 
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Consistently, the reconstruction procedure itself does not de- 
pend on the undetermined phase a(£, <j>) of the displacement 
operator, while the Wigner function does have such a depen- 
dence. The phase factor acts as a metric coefficient for the 
mapping from Hilbert space onto the phase space. We recall 
that for the harmonic oscillator, one can recover the Wigner 
function via an inverse Radon transform from the quadrature 
probability distribution 12711 . Equation d2"0t is then the analo- 
gous for our system. 

As a rather simple yet illustrative example, let us note that 
for the vortex state |^o) Ut is diagonal, so the tomograms 
co(4>', </>'/£) are independent of cf> and i and all of them equal 
to 1/(2-71"). Performing the integration in Eq. (f2Qb we obtain 
precisely the Wigner function (0. 

Finally, for the feasibility of the proposed scheme, we need 
the projection onto the eigenstates \<p). Since these states cor- 
respond to the measurement of a continuous variable, such a 
measurement can be done only approximately. A good ap- 
proximation seems to be the projection onto the wedge states, 
although other experimental schemes are also available 12811 . 
According to Eq. ( TT8l , this projection, together with the deter- 
mination of t, are the two main sources of measurement error. 
Afterwards, the reconstruction errors for each particular ex- 
perimental implementation should be evaluated following the 
standard recipes developed in, e.g., Ref. l29ll . 

In summary, we have carried out a full program for the 
reconstruction of generic vortex states, including a complete 
phase-space description in terms a bona fide Wigner function. 
Though the implementation of this scheme may differ depend- 
ing on the system under consideration, our formulation pro- 
vides a common theoretical framework on the Hilbert space 
generated by the action of angle and angular momentum. Ex- 
perimental demonstrations of the method in terms of optical 
beams are presently underway in our laboratory and will be 
reported elsewhere. 
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